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Abstract. The dynamics of a class of nonsymmetric gravitational theories is presented in Hamiltonian form. The 
derivation begins with the first-order action, treating the generalized connection coefficients as the canonical coor- 
dinates and the densitised components of the inverse of the fundamental tensor as conjugate momenta. The phase 
space of the symmetric sector is enlarged compared to the conventional treatments of General Relativity (GR) by 
a canonical pair that represents the metric density and its conjugate, removable by imposing strongly an associated 
pair of second class constraints and introducing Dirac brackets. The lapse and shift functions remain undetermined 
Lagrange multipliers that enforce the diffeomorphism constraints in the standard form of the NGT Hamiltonian. 
' Thus the dimension of the physical constraint surface in the symmetric sector is not enlarged over that of GR. In 

the antisymmetric sector, all six components of the fundamental tensor contribute conjugate pairs for the massive 
theory, and the absence of additional constraints gives six configuration space degrees of freedom per spacetime point 
■ in the antisymmetric sector. For the original NGT action (or, equivalently, Einstein's Unified Field Theory), the U(l) 

invariance of the action is shown to remove one of these antisymmetric sector conjugate pairs through an additional 
first class constraint, leaving five degrees of freedom. The restriction of the dynamics to GR configurations is con- 
, sidered, as well as the form of the surface terms that arise from the variation of the Hamiltonian. In the resulting 

^ ■ Hamiltonian system for the massive theory, singular behavior is found in the relations that determine some of the 

QQ ' Lagrange multipliers near GR and certain NGT spacetimes. What this implies about the dynamics of the theory is 

' not clearly understood at this time. 
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Introduction 



^ , The Hamiltonian formulation of GR (viewing the metric and the extrinsic curvature of a spacelike hypersurface 
O^' evolving in time as a first order system) has become popular in GR in recent years primarily due to various attempts to 
, canonically quantise GR. However the importance of the Hamiltonian picture goes beyond this formal procedure. In 
any classical field theory, the Hamiltonian formulation clarifies many dynamical issues of the system. The Bergmann- 
Dirac canonical analysis begins with a classical action, tells one whether it is consistent, and identifies the constraints 
L J and allowed choices of gauge. This identifies the physical constraint surface in phase space, the dimension of which 
rS gives the number of dynamical degrees of freedom in the system under consideration. The resulting first-order system 
(Hamilton's equations) is then ideally suited for numerical investigations of the initial value problem. As little is 
known of a general nature about solutions of the field equations in NGT and phenomenology tends to rely heavily on 
the properties of the static spherically symmetric solution^, it is hoped that some results of a more general nature 
will be accessible through Hamiltonian methods, or at the very least numerical investigations. 

In the case at hand, one is looking for a description of a class of Nonsymmetric Gravitational Theories (NGT 
will refer either to this type of theory in general or to the 'old', or massless theory, whereas mNGT will refer to 
the massive theory recently introduced in js^, |j) that will tell one which fields in the fundamental tensor are 
involved in the dynamical evolution of the system, and which (if any) are given by constraints on the initial Cauchy 
surface or left as freely specifiable Lagrange multipliers. The analysis performed here will treat the action for NGT in 
first-order form, treating the symmetric and antisymmetric connection coefficients as canonical coordinates, the com- 
ponents of the spatial part of the fundamental tensor showing up as (weakly equivalent to) the conjugate momenta. 
Hamilton's equations for the system are then similar to those of GR given in B] , with phase space suitably enlarged 
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^ |46[ gives the general form of the spherically svmmetric solution within the context of Einstein's Unified Field Theory, a special case 
of which was popular in the early works of NGT jsif ; another is the Wyman solution | ^^ , recently resurrected for both the 'old' |^ and 
massive [pi versions of NGT 
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to accommodate the antisymmetric sector, and the algebraic compatibiUty conditions left unimposed. Although it 
is possible in principle to treat the action as a second-order system (depending only on the fundamental tensor and 
its derivatives) by generalizing the inversion formula of Tonnelat p5| , p[ , the algebraic complications inherent in the 
solution of the generalized compatibility conditions make this approach unfeasible. 

The action for mNGT p7| , |3^ , Q has recently been introduced in order to make the antisymmetric sector well- 
behaved when considered as a perturbation of a GR background. The necessity of making this shift away from 
the original formulation of NGT was hinted at in some early works in NGT | p9| , p2| , but made a necessity by the 
work of Damour, Deser and McCarthy |^ in showing the bad asymptotic behavior of NGT perturbations about 
asymptotically flat GR backgrounds (reviewed in [Q). The action for niNGT is designed so that the skew perturbation 
equations have the form of massive Kalb-Ramond theory on a GR background, with additional curvature coupled 
potential terms. Thus in the case where the GR background is fixed, one finds three propagating degrees of freedom in 
the skew sector, and three algebraic relations that couple the remaining modes locally to the source (and background 
curvature), effectively removing what would be negative energy modes had they propagated. The results found here 
will indicate that in general, all six antisymmetric components of the fundamental tensor exist as independent degrees 
of freedom in the theory. 

The bulk of this paper (Sections §^) will be devoted to developing the formalism and recasting the dynamics of 
NGT in Hamiltonian form. This development is not particularly elegant, as the Hamiltonian is burdened by a large 
number of Lagrange multipliers, and the familiar patterns of the GR analysis become obscured. Where possible, 
parallels to the analogous GR analysis have been made, in order to not lose sight of the (in some ways very simple) 
overall structure. Indeed, although the algebra is sometimes rather messy, the analysis presented here is a relatively 
straightforward application of the Bergmann-Dirac constraint algorithm. 

Upon concluding the constraint analysis, in Section |^ the number of configuration space degrees of freedom per 
spacetime point in the old version of NGT is shown rigorously to be five, clearing up some confusion that has occurred 
in the literature ||2^. This however is not a new result, as it has long been known in the context of Einstein's Unified 
Field Theory . (Their analysis is not altered by the reinterpretation of the skew sector or the introduction of the 
source coupling into the action that defines NGT |Q.) The Cauchy analysis performed in the aforementioned works 
also determined that the Unified Field Theory (hence also NGT) 'suffers' from multiple light cones in that there 
are three causal boundaries within which different modes of the gravitational system propagate. It is expected that 
these 'causal metrics' will be slightly different for massive NGT, and in fact the considerations of Section || would 
seem to indicate that at least one of these has no pure GR limit (in the sense that as g[AB] 0, the related causal 
metric also vanishes) . The presence of more than one physical light cone confuses the issue of what constitutes the 
proper description of a Cauchy surface, and is discussed in more detail in Section 2.1. 



1. The Action and Field Equations 

The formalism developed in Section 3 of allows one to easily write the action for theories with a non-symmetric 
'metric' and connection coefficients in a surface adapted coordinate system suitable for the Hamiltonian analysis 
presented here. All covariant derivatives and curvatures have been defined with respect to a torsion-free connection 
(which in general is not required to be compatible with any tensor in the theory), whereas what is normally 
treated as the antisymmetric part of the connection coefficients is considered as an additional tensor A^^j. Despite 
the fact that this split introduces more objects into the formalism, the reduction to GR is made slightly more 
transparent, and allows the antisymmetric part of the torsion- free connection to be identified with the structure 
constants of a general basis in the standard wayQ. Although the connection will not initially be assumed torsion-free 
in the action, a tensor of Lagrange multipliers L^^ will be introduced to ensure vanishing torsion at the level of 
the field equations. In fact this is not necessary, as one could impose the torsion-free conditions from the outset 
and remove the appropriate connection coefficients from the action. This in fact is what will be done for most of 
the torsion-free conditions, as a judicious choice of arbitrary torsion terms added into the action causes the related 
connection components to disappear altogether. However, there are three conditions that have not been impos ed in 
this way, as it was felt that the formalism was simplified slightly by not doing so (discussed further in Section ^!l] ). 



^If one had instead defined a covariant derivative that was not torsion-free as in|g^, one would have to introduce additional fields 
into the action, and the split between the efli'ects of torsion and those due to the choice of coordinates would be obscured. 
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The action given in [Q {IGttG = c = 1) extended by the above mentioned Lagrange multipher term is then^: 
Each contribution to this will be denoted by corresponding Lagrangians (in order of appearance) , each to be separately 



decomposed in Section 3.1 : C^n,C^w, C^w^ : ^i, ^m, ^^t- Densities with respect to the determinant of the fundamental 



tensor g^s, will be indicated by boldface (eg. g'^^ := — det[g]g"^'^). The action (O) encompasses that for mNGT 
when a = 3/4, 'old' NGT (a = 0,m = 0) (equivalent to Einstein's Unified Field Theory), and recovers GR in the 
limit that all antisymmetric components of the fundamental tensor are set to zero. For simplicity, and to avoid 
the ambiguities inherent in defining the antisymmetric contributions to the matter stress-energy tensor, source-free 
spacetimes will be dealt with exclusively in this work. 

The Ricci-like tensor in the action has been split up into two contributions: R^g = Rab + Rab- "The first is 
identified as the Ricci tensor (i.e. it reduces to the GR Ricci tensor in the limit of vanishing antisymmetric sector): 

Rab = ecPs^] - esp^^] - ieA[rgp] -I- ^eBiT^c] + ^ba^ce - ^cb^ea + ^fAB]^EC^ (1.2a) 

and the second contains contributions from the antisymmetric tensor field A^^^j {A a '■= A^^): 

R% = [A]'Ab + Ve,, [A],, + AS^Agc- (l-2b) 
The Ricci tensor is defined from the two independent contractions of the geometric curvature tensor derived in the 
standard way from the connection that defines parallel transport of the general basis vectors Q by: Ve^ [ejs = F^^ec. 
The basis vectors are related in the usual way to a coordinate basis via a vierbein ca = Ea^B^, although the vierbein 
in this case is not necessarily mapping the coordinate basis into a Lorentz frame. As mentioned above, the torsion 
tensor: 

T^c = [Ve« [e]c - Vee [e]B - [es, ec]] = F^^ - F^^ - Cbc^, (1.3) 

will be constrained to vanish by variation of the Lagrange multipliers , allowing one to determine the antisym- 
metric part of the connection from the vierbeins. 

Variation of both the connection coefficients (F^^) and the antisymmetric tensor [A^q) will result in the com- 
patibility conditions: 

rrAB Y7 \„-\AB „DB \A „AD \ B , 2 x[A x^] „(DE)ij/ n /i a\ 

where equations (4.14) and (4.28) of Q have been written in terms of the densitised components of the inverse of the 
fundamental tensor. Note that the last term is precisely what is necessary in order for the compatibility conditions 
to be consistent with A a = 0. The remaining field equations are^ 

Aa = 0, (1.5a) 

1^ = iag^^^^VFs, (1.5b) 

Ve.[g][^^l=ag(^^>M^B, (1.5c) 

TZab = Rab + Ve,^ [W]^^ - ^aWAWe - {m'MAB = 0, (1.5d) 

where the mass tensor appearing in the last of these is: 

Mab = g[AB] - gCAgBDg^^"^ + 5gSAg'^^lg[CD]- (1-6) 



Neither (1.4) nor any of ( |1.5| ) will be used in order to pass to the second-order form of the action; they will instead 
be used in surface decomposed form to check that Hamilton's equations (as derived from the NGT Hamiltonian 
developed in the following section) are equivalent to the Euler-Lagrange equations quoted here. 

To date there has been no complete physical interpretation of the antisymmetric structure, nor will anything of 
the sort be attempted here. The action (^]|) will be considered as the fundamental starting point of the Hamiltonian 
analysis, and the goal of the present line of research is be to make some headway towards understanding the dynamics 
of the system. Thus in general one is left with many ambiguities as to how to make measurements on the skew 
sector, and the meaning of the various additional antisymmetric tensors as well as the additional contributions to 
the symmetric sector (GR) quantities is unclear. However, once one understands what data is configurable (Cauchy 



^For the translatio n of t he con venti ons used here for a coordinate basis with that of js^ see Section 1 of jij. 

*Note tha t one of (1.5a) and (|l.5q) is redundant, as has been replaced in the compatibility conditions derived directly from the 



action using (1.5h), to give the form (11.4). 
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data), the field equations fix the evolution uniquely (up to diffeomorphisms), allowing one to determine the physical 
implications of the presence of the antisymmetric sector in cases where the interpretation of the symmetric sector is 
essentially that of GR. This will be the case in the asymptotic region of spacetimes that are asymptotically dominated 
by the symmetric sector, as well as in perturbative scenarios (discussed in Section |^) where initial data is close to 
GR configurations. 

The next step in the analysis is to introduce the foliation of the spacetime manifold M, and decompose all objects 
into components perpendicular and parallel to the surfaces. 

2. Metric and Compatibility in a Surface Compatible Basis 

Normally one would assume at this point that spacetime is globally hyperbolic, so that a Cauchy surface Eq 
exists in M (a closed achronal set without edge, intersected by any smooth causal curve that is future and past 
inextensible). Spacetime may then be viewed as the evolution of 3-geometries, that is, (M,g) emerges from the set 
of field configurations on spacelike hypersurfaces that make up a foliation of spacetime: {(St, gt)\t G I C M}}. Then 
one may show that diffeomorphically equivalent initial data on Eq generate spacetimes that are diffeomorphically 
equivalent, and a physical spacetime is then identified as an equivalence class of such solutions (equivalent up to 
diffeomorphisms) . 

For NGT however, a sensible definition of hyperbolicity is somewhat more complicated. In Chapter IV of [^, it is 
demonstrated that the field equations of Einstein's Unified Field Theory propagate information along three separate 
light cones (this also occurs when considering covariant wave equations for particles of spin greater than or equal to 
one 10). These cones are defined by symmetric metric tensors, the set of which will be denoted: {gc}, and referred 
to as the causal metrics of NGT^. Thus the statement that a vector is timelike (spacelike) with respect to {gc} will 
indicate that it is timelike (spacelike) with respect to all of the causal metrics in {gc}. Global hyperbolicity would 
then require that Sq be achronal with respect to {gc}, and the causal curves would of course also be causal with 
respect to {gc}- This would ensure that none of {gc} would become degenerate anywhere in spacetime, and hence 
that So is in fact a Cauchy surface. At this stage it is not known whether one can require such a condition on the 
fundamental tensor^. However even if this is not possible, the Hamiltonian system is still relevant for considering 
the dynamics of the theory locally. One may instead specify initial data on a large enough 5*0 C Eg in order that the 
configuration on St C Ej (St is a closed achronal set) is determined uniquely. This requires that {gc} is known in 
order to determine i)+(5o) or D^{St) (the future (past) Cauchy development or domain of dependance of So (St)), 
to ensure that So C D+(S'o) or So ^ D^{St). For the purposes of this work, it will suffice to discuss the global 
problem, as one may restrict the results to in a straightforward manner once the causal metrics of mNGT are known. 

2.1. The Surface Adapted Basis. We begin by assuming that there is a time function t that is used to foliate M 
into hypersurfaces of constant time Et. This requires that each (3-)surface of constant time, Et, be spacelike with 
respect to {gc}, and Vt* = 1 defines a vector t that is timelike with respect to {gc}. This ensures that the degrees 
of freedom that travel along each of the light cones in {gc} will evolve forward in what has been chosen as time. 

In GR, one may then introduce a coordinate basis on the surface and a unit normal, such that the metric is of 
the form: g = 0-^ 9-^ — ^abO'^ ® 0^ . This effectively reduces the metric to a Riemannian metric on E, and the 
projection of dt perpendicular and parallel to the surface are the lapse and shift functions respectively [|l6|. It is the 
components of this metric on E that are taken to be the Canonical coordinates in the Hamiltonian approach to GR, 
and all other spacetime tensors are decomposed onto the surface by considering the components perpendicular and 
parallel to the surface as separate tensors on S. 

On attempting to generalize this to NGT, one finds that the presence of more than one physical metric on M 
implies that there is no physically well-motivated or natural choice of metric on E that would play the same role. 
One may introduce a coordinate basis on the surface in the usual way, but in general each of the metrics in {gc} will 



^In the Unified Field Theory, {gc} = 7} where I represents gt^*") and its inverse, h represents g(^^) and its inverse, and 7 and 
its inverse are defined by 7"'' = i^h'^l^ — I'^P where h = det[hp,y] and g = det[g^i,]. These metrics were found to be compatible in the 
sense that there is a largest and smallest speed of light, and all three metrics merge into the Riemannian metric of GR in the limit of 
vanishing antisymmetric sector. As the set of causal metrics is not known explicitly for mNGT at this time, strictly speaking one cannot 
discuss the causal structure of any given spacetime. 

^What is commonly referred to as the 'metric' of NGT is the tensor determined directly by the field equations (appearing in (1.1)), 



the determinant of which defines the spacetime volume element in the action, and from which one derives the causal metrics {gc}- The 
presence of antisymmetric components indicates that it is in fact not a metric (in the usual sense) at all, and will therefore be referred 
to as the 'fundamental tensor' (in accordance with Lichnerowicz p7|). Note however that one may consider it as a Hermitian metric in 
a complex or hyperbolic complex space j^, ^sj , or in more general spaces ps) . 
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provide a different definition of the unit normal vector e±, and the resulting decomposition of tensors onto S will 
be inequivalent. However as this is merely a choice of parameterization, the physical content of the system will be 
independent of the choice of metric used to define the surface decomposition. In this work, the components of g'^^^) 
will be used to define the unit normal e±, therefore making a particular choice of decomposition: 

g-i = e_L ® e_L + B^ejL d) Ca - B^Ca ® ej_ - -f'^^Ca ® et, (2.1a) 

gsy'mm = 6^ ® 6^ - T^'^^^Ba (g) Cfc, (2.1b) 

where 'symm' indicates the symmetric part of the inverse of the fundamental tensor. This choice is made in order 
to simplify as much as possible the form of the action (1.1), as the fundamental tensor appears in the form g"*^. 

The time vector is decomposed as t = t^BA = Ne± + N°'ea, where {N, N"") are the lapse (component of t along 
e^) and shift (projection of t on the surface Sj). This basis can therefore be written in terms of the coordinate basis 
{dt,da) by: e± = j^dt — ^da, = da, and the covector bases by: 9^ = Ndt, 9°- = dx"- + N°-dt. This defines the 
vierbein through 9^ = ^dx^, the determinant of which is: 

E ■.= dei[E'^^,]= N. (2.2) 

Note that is a coordinate basis on S and can be therefore be thought of as acting on tensor fields as an ordinary 
derivative da- It is also useful to denote surface vectors (not just the components) by V := V^Ca and covectors by 
Q :— uJa9°' . It is not difficult to check that in the coordinate basis {dt,da), g^mm ^^'^ its inverse take on the usual 
ADM form g. 

The fundamental tensor is found from the inverse of ( 2.1a ): 

g = F9^<E>9^- {aa - I3a)9^ ® r - (a, + /3,)r ® 9^ - Gab9'^ ® 9\ (2.3a) 

where: 

F := 1 + G(^ab)B°'B'' , aa := G[afc]-B^, Pa ■= G(^ab)B^, (2.3b) 
from which one finds UaB'^ = 0. The spatial part of the fundamental tensor {Gab) is given as the inverse oi -f — B (g) B: 

_ B''B'^)Gct = GbciY" - B'B'^) = 5t. (2.3c) 



Note that invertibility of the fundamental tensor ( |2.3D is required in order that the volume element in the action be 
nondegenerate, and this requires that Gab exist. Given that the fundamental tensor is a nondegenerate solution of 
the field equations in no way guarantees that {gc} are all nondegenerate, allowing for the possibility that a regular 
solution of the field equations has regions in spacetime where NGT perturbations 'see' a singular spacetime. 
Some useful identities may be derived from ( |2.3c| ): 

Varying ^/—g with respect to the densitised components of the inverse of the fundamental tensor results in: 

= -j^PaSB^ + ^GbaSr'- (2.5) 

This will be used order to compute the variation of the Hamiltonian of NGT with respect to the same densitised 
components, as well as applied as an ordinary derivative relation in the surface-decomposed compatibility conditions 
in Appendix A.l. Note that the condition F = 2 corresponds to the case where \/—g is independent of (B'',7"''), 
and F ^ 2 will be assumed throughout this work. 

The non- vanishing structure constants for this basis are found to be: 

[e±,ea] = Cia^e^ + G_i_a''eb, 

C^a^ = CaMN)], C^a' = j^ea[N% ^^'^^ 



and if one were to require that the connection T be torsion- free using (1.3), the antisymmetric components of the 
connection would be: 

^ta] = lea[HN)], Tl^^^A_ea[N\ rf,,]=0, rf,,j=0. (2.7) 
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A final useful property of this basis is: 



0, 



(2.8) 



which can either be derived directly from a contraction of the Jacobi identity, or proved by a brute force insertion of 



the structure constants (2.6). This will b e us eful when identifying the diffeomorphism constraints, linking them to 
the algebraic field equations in Appendix A. 2. 



In beginning with the components of the inverse of the fundamental tensor in (2J^), the configuration of the system 
is being described by the degrees of freedom iV, A^"^, B°,7°^ (the latter two are included as densities since they will 
turn up later on as weakly equal to the conjugate momenta in this form). As one can see from the form of ( ^.3| ), 
different parameterizations of the fundamental tensor will have a large effect on the details of the analysis presented 
here. If one had begun by requiring that the symmetric part of the fundamental tensor took on the ADM form, 
then the inverse (that appears in the action) would have had nontrivial g^^") components, and the analysis would 
have been more complicated. In NGT, one must be slightly more careful to be consistent in what one is calling the 
independent degrees of freedom, as the antisymmetric degrees of freedom will mix with the symmetric sector as one 
moves from spatial components of the fundamental tensor to the spatial components of its inverse. 



2.2. Surface Decomposition. Four-dimensional covariant derivatives are decomposed as usual mM into surface 
covariant derivatives (written as Va ) whose action on basis vectors in TS is given by Va [e]f, — r^^Cc, and 
contributions from derivatives off of E which will be written in terms of the surface tensors: 





r 


= r 






= r 


^ ±6j 


< 


= r 



a 

kab := r, 



ra 



(2.9) 



ab- 



(Both k and F^j, will initially be considered as nonsymmetric, allowing for the possibility of non-zero torsion. The 
torsion- free conditions will be imposed for these components very early in the development, and from then on 
they will be assumed to be symmetric.) In GR most of these are related by algebraic compatibility conditions 
(cr" = ^"■^ah.Ca = 0,F = 0,u\ = ^"''^kbc), and the absence of torsion would allow one to further relate a and {u,w) 
to the structure constants. The skew tensor A: 



the torsion Lagrange multipliers: 



A 



lv 



]ab 



L 



ah 



A 



L 



ah 



A° 



r a J 
^b 



^bc 



AO 



'- be 



L 



be 



(2.10) 



(2.11) 



and the vector fields: Wa = (VF, Wa) and V 



w'^, and the traceless part: 



Tb 



- {IJ""), are all decomposed in a similar fashion. Traces are denoted as 
l/3S^v. The components of A a are given by A^ — v,Aa — ba + 



6' 



and jab is an antisymmetric surface tensor by definition. It also will be useful to define a few combined quantities 
that will appear: 



Kab 
K J, 

J b 



V 



— kab + jabi 

= hil'^Kte + 7'°i^cb) = I'-^'hbe + -/^"""^Jbe 



(2.12a) 
(2.12b) 
(2.12c) 
(2.12d) 
(2.12e) 



Aside from Kab, which will be useful as a notational tool to combine the results of both sectors into one relation, 
these combinations correspond to symmetric (kab, u^) and antisymmetric {jab, v^) sector contributions once an index 
has been 'raised' by the spatial fundamental tensor 7°^. 



2.3. Lie Derivatives. The standard Lie derivatives defined on E: 

£%\y] = [A,y] = {X-ea[Y'']-Y-ea[X%e,, £f[Y^]=X'^ea[Y^i 



X 



£f[Q] - (A"eaM +c^aee[A''])r, £'^'[^A - ^"e^K], 



.(3)r 



(2.13a) 
(2.13b) 
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(and the standard generalization to higher-order tensors) will be used throughout. Note that the perpendicular 
components are treated as scalar fields as far as surface defined derivatives are concerned. This is also true of the 



derivative off of S 1 16 



Ax^eAy]^X^{e^[Y'']+Y''C^a^)eb, d^^eAY^ = X^e^Y^ (2.14a) 

dx^eM = X^(ei[c^6] +co„Q/)^?^ dx^e^K] - X^e^K], (2.14b) 

where the structure constant takes into account the possibility that the surface basis may 'move' on S as one moves 
perpendicularly off the surface. These definitions are also extended to arbitrary tensors in the standard manner. 
From the spacetime Lie derivatives of tensor densities: 

£x[T] = ^/^£x[T]+T£x[^/^. (2.15a) 

£x[V^ = ig^^i^xfesd = -\^CB£x[g''% (2.15b) 
and using the defined surface Lie derivatives, one finds: 

dt[A^ = dt[^/^ := iVdeJx/^ + £f[A^ = N{eAV^ - T^Cla") + Vi^) [N]». (2.15c) 

Time derivatives of arbitrary tensors or tensor densities are given as in GR by: dt[Z] = dt[Z] = Nd^^ [Z] + £^i^ [Z], 
where Z is any tensor density. 

One would now like to identify the extrinsic curvature of the surface S under consideration. Although there are 
a variety of equivalent ways of describing it in the context of Riemannian geometry, none of these are equivalent 
for the type of nonsymmetric theory considered here. One way is to compute the perpendicular component of the 
parallel transport of a surface vector along the surface (Gauss' formula 0): 0^[Vjg[N]] = T^^M'^iV'' = katM^'N'', 
in conjunction with the fact that the surface projection is just the surface transport result in GR 0°[Vjg[7V]] = 
6'" [V^^' [iV]] 1^. One could equivalently describe it by the change in the unit normal as it is parallel transported 
along the surface (Weingarten equation) Va[e^] |Q. In GR, the magnitude of this vector does not change when 
transported along the surface (since T-^^ = Ca = 0), and the result is a vector tangent to the surface describing the 
amount that the normal vector must 'rotate' in order to remain normal VQ[e_L] — ^a±^b — u^a^b- Metric compatibility 
relates these two definitions since kab — Jbcu'^a- the case of NGT, this definition is clouded by the fact that in 
general Ca ^ 0, and the length of the unit normal will not be preserved under (this definition of) parallel transport. 
One could also identify the extrinsic curvature of S with the change of the 3-metric perpendicularly off the surface 
Ve^ [g][ea eb] = -de^ [Aab + lcbU% + lacU% = 0' giving de^ [Aab = 2fcab (or the equivalent derived from the inverse 
metric) . 

Although equivalent in GR, none of these definitions coincide in nonsymmetric theories in general. Indeed, since 
the compatibility conditions cannot entirely be written in terms of a metric compatible connection, there is no 
truly natural choice of parallel transport, and any choice of 'extrinsic curvature' based upon a covariant derivative 
is conventional at best. The approach followed here will be the first of those described above, essentially due to 
convenience. It turns out that with the choice of decomposition made here, F^j, — kab is a canonical coordinate, and 
as such plays a more central role than u'J, (which will turn out to be a derived quantity, as it acts as a determined 
Lagrange multiplier), or the relation derived from the evolution of the spatial metric. 

The stage is now set, as one now has the tools necessary to decompose spacetime tensors and covariant derivatives 
into spatially covariant objects, and identify time deriv ativ es in the action. The configuration space variables have 
also been chosen to be: [N, iV"°, B°, 7°^) . In Appendix |A.l| the compatibility conditions (F4) are decomposed, and 



the field equations (^]^) in Appendix A. 2 . In both cases, algebraic conditions that will determine Lagrange multipliers 



on E have been identified, as well as those responsible for time evolution. The reader is again reminded that all of 
these will be derived from the Hamiltonian, and are only given here in order to introduce some defined quantities 
(£,25,3) that will help to simplify the presentation of the algebra, and ensure that one is finding results that are in 
accord with the Lagrangian variational principle. 



3. Determining the Hamiltonian 

This section is concerned with obtaining the form of the Hamiltonian for NGT. The various terms in the Lagrangian 
( |l . l[ ) will be decomposed into surface compatible form, and the fields that appear as time derivatives identified. After 
this, writing down the form of the Hamiltonian turns out to be a rather simple task, however the work is far from 
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done. The Hamiltonian is riddled with Lagrange multiphers, and the constraint analysis is where one begins to see 
the full structure of NGT emerge. 



3.1. The Surface Decomposed Action. The simplest Lagrangian density terms in (1.1) are decomposed as: 



Ci=N\v^^ + Nr{ba + \l,), 

C^W = A^B° [Vi^) [W] - de, [W]a + W{aa -Ca)\+ N'^^-''^Vf^ . 



(3.1a) 
(3.1b) 
(3.1c) 
(3.1d) 
(3.1e) 



The first of these defines Tim, which will be useful in Section 4.3. The remaining contribution /Ijjns, will be further 
split into contributions from each of the symmetric and antisymmetric sectors as and Cj^a, and further split 
to give contributions from each component of the tensor (for example, >C_R(„j) corresponds to the contribution from 
N ^'^"■^^ Ri^ab)) ■ Note that the fact that g(-'-°' = in this decomposition implies that there is no contribution from 
^R(±a) ' which is essentially why the choice of surface adapted basis (2.1) was made. Making a choice of parameteri- 
zation in which g^-"-") ^ would have explicitly introduced a term involving dtB°' into the action (see (A. 5b)). The 
remaining contributions are: 



•■[ab] 



c 



iVl^''') [bah ~ JacV% - JbcVl + A^Ag,] , 

-iVB'^ [de, [A]^, - de, [b]a + 2Vf ^ [v]\ - V^^^) [v] 

+ 2bbu\ - 2bau + 2(at - Cb)v\ - (a^ - Ca)v - 2u\Xl^ + 2jab(J^] , 

- jacU^b - jcbu'^a + (T + u)jab - KaV^ + kcbV^^ + aabb - ttbba + acX^h] , 

-NV^[^de^ [u] + Vi3)[a]- + Tu + a''{aa - 2ca) - ulu"^] , 

^7^°'' [R% + de^ [%ab) - [a]b) - daab + (r + U)kab ~ k.bu'a - kacU%] , 



(3.1f) 



where R(^Ji,-^ is the surface Ricci tensor given by (1.2a) determined from surface covariant derivatives, and the identity 
R[AB] = (see Section IV in |Q]) has been used. 

The terms that result in time derivatives in the action are then: 

^NB-d,^ [W]a - NB-de^ [Xtb + NB-d,^ [b]a + NV^de^ [u] + N-y'^'de^ [KU, 

and so defining the vector field: 

Wa^-Wa-Xlb + ba, 

will simplify the remaining analysis. Rewriting the above in terms of the time derivative gives: 



(3.2) 



-B''£f[W]a 



.,ab r 



where: 



£f [K]ab = iV^V(3) [K]ab + i^cbVi^) [N]- + KacVf^ [N] 



(3.3a) 
(3.3b) 
(3.3c) 



Although these have been written in surface covariant form, all dependence on the spatial connection coefficients 
vanish in these. As it is clear that W a will be a canonical coordinate, (3.2) will henceforth be used to replace Wa 
everywhere in the remainder of this work. 
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Gathering together terms in the action that are multiphed by the same component of the fundamental tensor and 
identifying the time derivatives as above, the Lagrangian density for NGT can now be written: 

+ \a{Wf ~ vlv\ - V(3) [af ~Tu^ a^iaa - 2c,) + u%u\] 
+ NB'^ [V(3) [W] + Wiaa -ca)- 2Vf ^ H", + V'^) [v] 

- 2hu''a + ^bau - 2(afc - Cb)v'', + (fla ' Ca)v + 2u\\l^ - 2j,ba''] (3.4) 

- iacU% - jcb" a + (r + U)jab - KaV^ + KhV^a + - O^bba + O-cKb] 

+ iYit-^^) [RfX) - Vf^ [a]b) - a.flb + (F + u)Kb - fccb"^ - ^ac^^ 

- iaC^a + A'^,, - baWb + Atd - fefc) + bah ^ JacV\ - JbcV'a + KcKd] ■ 

Up to this point, the torsion- free conditions have been apphed arbitrarily in order to simplify the Lagrangian density 
as much as possible. At this stage, the fields k[ab] and F^^^j have been dropped along with their associated torsion 
Lagrange multipliers L"** and L"**, as they would appear only in these torsion Lagrange multiplier terms. From 
this point onwards, both kab and T^f, will refer to symmetric objects. The Lagrange multiplier term imposing the 
torsion-free condition: 

w% = u% + Clb", (3.5) 
has also been dropped, as once all terms (in the action as well as compatibility conditions and field equations) are 



written in terms of dej^, will not appear either, and may be determined by (3.5) once u'\ is known. One would 



also like to impose the torsion conditions to solve for explicitly, however this would introduce a number of terms 
involving the structure constants C±a^ into the action, compatibility conditions and field equations, and it was felt 
that avoiding this simplified the algebra. Although there is nothing preventing one from doing this, the presentation 
is slightly simpler if the torsion condition is not imposed in this case. 



3.2. Legendre Transform. In the form (3^), the fields with conjugate momenta are trivial to identify (these are 
densities by definition, and will not be written in boldface)^]: 



S[dtu] 

5[dtWa] 



-g, (3.6a) 
B^ (3.6b) 



'^'^"^=#^1-^"'' (3.6c) 

0[dtKab\ 

where the action has been written in terms of the Lagrangian as: S'ngt = J dtLfiCT- One could in principle treat 
all fields in the action as canonical coordinates, generating many cyclic momenta and hence many second class 
constraints that are trivially removed by the introduction of Dirac brackets Instead the more economical 

approach of treating any field in the action whose time derivative does not appear anywhere as a Lagrange multiplier 
field will be followed^. 

The extended Hamiltonian density is written as: 

n,=no + Ncj^abiTT'"'' - r') + 2N^a{p'' - B°) + N d){p - (3.7a) 



where the form of the Lagrange multiplier terms that enforce the algebraic conditions (3^) have been written as the 
reduction of the covariant form: E(f>AB{T^^^ — g"*^)- The basic Hamiltonian density is derived from the Lagrangian 



^Note that the Canonical coordinates have been taken to be connection components (and Wa), resulting in the densitised components 
of the fundamental tensor appearing as conjugate momenta. By discarding a total time derivative in the action, one could easily arrange 
for the fundamental tensor degrees of freedom to play their more traditional role as Canonical coordinates. 

*It is common to include the lapse and shift functions and their conjugate momenta in phase space, particularly if one is attempting 
to realize all of the generators of the spacetime diffeomorphism group on phase space |l7| . This will not be done in this work simply 
to avoid the unnecessary algebra, as it is a trivial matter to extend the formalism later if necessary. 
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density in the standard manner: Tig = X)/ P^Qi ^ ^^igt, and is found to be (equivalent to dropping time derivatives 
in (3.4) and taking its negative): 

Ho ^Hrn + V^£f M + B^£f [W]a + l"" £f [K^ 

- N^[lv + l-{ba + + 2i"(Ca -aa + C^a^) 

+ h^iWf - vlv\ - Vi^) [af -Tu~ a'^iaa - 2c,) + u\u\] 

- iVB'^ [Vi^) [W] + W{aa - ca) ~ 2vf H^, + V^^) H 

- 2hu\ + 2baU - 2(06 " Cb)v\ + {aa - Ca)v - 2jab<J^ + 2u^,AgJ (3.7b) 

- jacU\ - jcbu''^ + (r + u)jab - kcaV% + kcbV"^ + Gabb - abba + a^Xlb] 

- 7V7('^'') [r'^^^^ - Vfj [a]b) - aaab + (T + u)kab - Kbu'a - kacU% 

- ^aiWa + A",, - ba){Wb + X'id - bb) + babb - jacV^b " JbcV^a + >^tKd] ■ 

Hamiltonians and Hamiltonian densities will be related everywhere by (for example): H^, = (PxTi^,, where H will 
refer to the Hamiltonian, and TL to the related Hamiltonian density. 

The canonical phase space variables that describe the configuration of the system at any time are: 

{{P',Ql)} = {{p,u),{p\Wa)A^''\Kab)], (3.8) 

and the usual canonical Poisson brackets have been assumed: 

{Qi{x),Qjm - {P\x),P-\y)} = 0, (3.9a) 

{Qi[f'],P'[gj]}= I d'zf\z)gi{z). (3.9b) 



The notation Qi[f^] is used to denote the smearing of the object by an appropriately weighted test tensor (eg. 
P°'[fa\ — d!^ zp°-{z)fa{z)) leaving a scalar on S, and the sub/superscript T indexes the sets in ( |3.8| ). The Poisson 
brackets are standard: 



{F[Q,PIG[Q,P]}= f d'zJ2 



6F SG 6F SG 



5Qi{z) SPi{z) SPi{z) SQiiz) 



(3.10) 



where {F, G) are scalar functions of the canonical coordinates. 

The fields (p, u), although symmetric sector fields, do not show up in the standard Hamiltonian formulations of GR 
as canonical variables. There one has u = ^""^kab, and this is used in the action to either remove the time derivative 
from u (second-order action) or write: dtu — ^""^dtkab + kabdtj"''' [0. The point being that in nonsymmetric 
theories, one cannot replace all terms in the action by various combinations of the spatial fundamental tensor and 
extrinsic curvature in a straightforward manner, and it seems one must be content to carry around the 103 Lagrange 
multipliers: 

{N}:^{N,N''}, {L}:={iL-),il,n}, m :^ {<P, ^a, ^(at), ^[ab]}, 
{7}:={B^7'^^}, {A}:^{W,ba,v^Xl}, {F} :- {F, c„ a,, a^ z«1„ F^J, (3.11) 
where the Lagrange multipliers have been split up into smaller subsets for later reference. 



4. Constraint Analysis 

Now begins the analysis of the constraints and conditions that follow from variations of the Hamiltonian with 
respect to the Lagrange multipliers. It will be found that all of the compatibility conditions and field equations 



(1.3 1.4,1.5) are reproduced from the Hamiltonian, leaving {N} as undetermined Lagrange multipliers enforcing the 



diffeomorphism constraints in the standard form of the Hamiltonian. Two second class constraints will also be found 
which could be used to strongly remove the conjugate pair {p, u) from the phase space of the theory, leaving the 
symmetric phase space with the same dimension as in GR. As the results are unsurprising and the algebra not unduly 
difficult but rather lengthy, the results will be given as sparsely as possible. 
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4.1. Constraints Enforced by {</>}, {L}, {A} and {F}. The Lagrange multipliers {</>} are designed to enforce the 
conjugate momenta conditions ( |3.6| ): 



5H^ 



2iV[p° -B"] K. 0, 



(4.1a) 
(4.1b) 
(4.1c) 



These conditions clearly determine {7}, with one constraint left over: tpi^ which has been identified in (4.1c). From 
{i}, the Lagrange multiplier enforcing the remaining torsion constraint from (2.7) gives: 



5Ih 
(5L 



^^-2N[ca-aa + C^a^] « 0, 



which will be used to determine Oa- The constraints resulting from the variation of are: 

5H^ 



5\ 



-Nv « 0, 



(4.2) 

(4.3a) 
(4.3b) 



These are the components of (1.5a), and will be used to determine v and A^^ respectively. 

The conditions arising from the variation of {A} are now calculated, first varying with respect to W and then 
„h ■ 



5H, 
5W 



Using the result that v ^ from ( 4.3a ), the trace of the last of these gives 1 ~ iVc'^''[B]'^ 
Next varying ba gives: 

^ « 27ve:f 1 -N\r- ^aj^^'^Wh] - o, 

Oba o L c, 

and using the above conditions one finds that: 



SH, 



(4.4a) 
(4.4b) 

^W, so that €1^"^ w 0. 

(4.4c) 

(4.4d) 



ab 



which, combined with the previous result, yields 1" « ^'y'- Wt,- At this point, all of the skew sector algebraic 
compatibility conditions ((A. 2c) and (A. 2d)) as well as the field equations (1.51) have been reproduced. 
Varying {F}, one finds: 



6H^ 
ST 

SH, 

SCa 



TVCi^") =: NiP2 « 0, 

-27V [L'^ + fj" + B''v1 - iB'^M/] w 0, 



(4.5a) 
(4.5b) 
(4.5c) 



whe re the first of these yields a constraint (which has been identified as '^2 in (4.5a)) discussed further in Section 
and the last is used to determine the Lagrange multiplier L°' as: 



4.2 



Using these, one computes: 



SH, 
Saa 



1 



b -r- ^B'^W. 



(ab) _ g-(-Lo)j 



(4.5d) 
(4.5e) 
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and the trace- free part of the variation with respect to u''^ gives the variation with respect to u^xa- 



(4.5f) 



which yields C^"*""^ « when the constraint (4.5a) is taken into account. Finally, the variation with respect to the 
surface connection coefhcients yields: 



6H, 

6r 



(4.5g) 



ab 



and all but €t of the symmetric sector algebraic compatibility conditions (( A.lc ), ( A.le ) and ( A. 11 )) have been 
found, and the Lagrange multiplier L" determined by ( d.Sd ). 

4.2. Time Evolution of the Canonical Fields and the Constraints Time evolution of the canonical 

variables is determined as usual (from (3.9) and ( 3.10| )) by Hamilton's equations: 



dtQi^{Qi,H,} 



dtP' = {P',H,}^- 



SIh 
SQi 



(4.6) 



These evolution equations should reproduce the same dynamics as given by (1.3, 1.4, |1.5| ). It is convenient at this 
point to show that the canonical momenta evolve in accordance with the dynamical compatibility conditions, and 
that the evolution of the canonical coordinates also allows one to identify conditions that must be found in order 
to reproduce the NGT field equations ( 1.5d| ). This is why this calculation is performed at this stage, even though 
strictly speaking it is not part of the constraint analysis. 

The functional derivatives of with respect to the canonical coordinates will be necessary: 

/) H 



5u 

SWa 
Skab 
6jab 



.^(-^b) V(3) [NY + 7^'^'=) V(3) [N]'' + 7^ 

-iV[7("'')(r + 'u) -2u"''], 
,^[afc]y(3) [iv]= +7['''=lv(3)[iV]'' +7[^''lv(3)[Ar]'' - iV^v(3)[7][H 



■iV[B"c 



BV 



bn 



)'["''i(r + u) + 2v' 



abl 



(4.7a) 

(4.7b) 
(4.7c) 
(4.7d) 



Time evolution of the conjugate momenta can be shown to be weakly equivalent to the dynamical compatibility 

5H, 



conditions (|A.l4 |A.l4 |A.24 |A.2bD: 



dtp = -- 



Su 
dtp = 



dt[7] 



ab 



ab 



(4.8a) 
(4.8b) 

(4.8c) 



since in all cases (heuristically) dtP ~ dt[7], and one is left with the result that €± ^ in evolution. Using ( [A.5| ) 
and the evolution of the coordinates: 

r rr 

dfu = -A « N(j), (4.9a) 



Sp 

a,w. ^ 



dtKab = 



Sp'^ 

§^ab 



' N<j)ab, 



(4.9b) 
(4.9c) 
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allows one to define the objects {91} ■ 



1^2 

4 



N 

.(3)r 



y^ab := N(bab - £%'[K]ab + N^ab " ^m^NMab- 



(4.10a) 
(4.10b) 
(4.10c) 



In order for Hamilton's equations (4^) to reproduce the dynamical field equations ( 1.5d ), one must find that {$K} ~ 0, 

ensuring that {0} are correctly determined. 

There are two constraints {■0} that have appeared from the variations so far (( 4.1c ) and ( 4.5a ) respectively): 

^l«p~y^[p^^'^^]«0, (4.11a) 
^2~pu-n''''Kab~0, (4.11b) 



where the density has been written as functionally dependent on the other momenta in order to stress that it must 
be considered as a functional of those fields alone (as in (2^)) when calculating Poisson brackets. Although all of 
the Lagrange multiplier variations have not been dealt with yet, the conditions that result from requiring that these 
constraints are preserved in time will prove useful at this stage. Requiring that ipi be preserved in time results in: 



5H^ 



2 SH, 



1 „ 



2~F'^''SWa 2-F 



2-F 



N€t « 0, 



(4.12) 



where ( [2.5[ ) has been used. This directly gives the remaining algebraic compatibility condition (A.4a). The preser- 
vation of "02 gives: 



5*02 = 



6u 



'SK^b 



F~2 



(4.13) 



where is defined in ( |4.10| ) and ©-^ by dA.lOa] ). 



4.3. Variations of {7} and {N}. In order to compute the variations with respect to {7}, it is convenient to rewrite 
the extended Hamiltonian density (3.7a) in the form ({3} is defined in (|A.^)): 



n,^nm + V^(iV3±± + i^g^N - iV0) + B''(-2A^3[a±] + £%'[W]a - 2iV0a) 



(3) ru/i 



+ r'i-N^ab + £f [K]ab - N^ab). (4.14) 

where ~ refers to the fact that some conditions and constraints have been imposed in 7i* , provided that the results 
of the variations are not weak ly aff ected. (Note that this form holds for variat ions with respect to {7} only.) The 
variation of TLm as defined in ( [3.1a| ) can be written (once again making use of (^.Sj)): 



l3aM±i_)5B'' + -Nm'{Mab- 



F 



GbaM^i_)S'Y'''' 



and one then computes the conditions enforced by {7} to be (using ( 4.10| )): 

5H, 2 



SB'' 

g^ab 



'^ab 



2-F 

1 



2-F 



/3a9^ « 0, 
Gba^ ~ 0. 



(4.15) 

(4.16a) 
(4.16b) 



As it stands, the Hamiltonian H^, is not in the standard form: j^{NTi + N"-Tia)- It may be put in this form by 
the removal of the surface term: 



Ft, 



dSaWWbB" + 2iV''k'';, - 2N1,'' 



as 



so that 



H* = j d^x{Nn + N^Ua) + F, 



(4.17) 



(4.18) 
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(Further discussion of surface terms will occur in the next section.) One then finds directly the Hamiltonian constraint 
v.: 

S TT 

n:=j^^ -20\ « 0, (4.19a) 

which has been weakly identified with the algebraic field equation ( |A.10a ) by making use of previous results. This 
constraint, combined with (4.13) and (4.16), establishes that {91} « 0. The momentum constraints are: 

ebip'Wa] + 7r'''ea[Ki,c] + pea[u] - 2eb[k\] 



gv(3)M-2vp)[k]'', 



(4.19b) 



0, 



which has been identified with the related algebraic field equations (A. 10b). The identity (A. 7) has been used in 
order to make the identification of Ti. and Tia with <&\ and respectively. The form of the Hamiltonian constraint 
is not given explicitly in (4.19a) as it can essentially be read off from (3.7a) or (A. 10a) directly. Unlike the case of 
the momentum constraint, there seems to be no relatively simple formulation of the Hamiltonian constraint in terms 
of canonical variables alone. This leads to computational problems (as will be discussed in Section 4^), as Poisson 
brackets with 7i cannot be calculated in a straightforward manner. 

With these last constraints, the Euler-Lagrange equations of NGT have been fully reproduced by Hamilton's 
equations combined with all of the algebraic compatibility conditions and field equations resulting from the variation 
of the Lagrange multipliers. 

4.4. Surface Terms. Up to this point all surface terms have been consistently ignored. This is fine if S is closed 
((9S] = 0) and there is no surface to have terms on, but in general one must add a surface contribution to the action 
(1.1) designed to cancel off any variations of i?* that do not vanish on the dY,. This ensures that the field equations 
that have been generated follow correctly from the Hamiltonian |40|] . 
To begin with, one may add to iJ, the surface contribution: 



JdY. (4.20) 
+ 7V7r['"=lA|J^ + 2NTT^''%b - 2Np^vl + Np^'v + Np^W], 

which ensures that the variation with respect to any of the determined Lagrange multipliers never result in any 
surface terms. In addition, there are surface terms involving variations of {7} and {N} to consider, resulting from 
variations of -Elm as well as contributions from the variations computed earlier in this section. These may be split 
up into terms involving the lapse and shift functions separately, and the full Hamiltonian for NGT including surface 
contributions written as: 



H^gt:= I dPx{N'H + N''Ha)+E,+E^M + EN + E^, 



where E^, is given by (4.17). The contribution from the shift function must satisfy: 



5E^ 



dSa[N''pSu + N''7r'"'dKb, 



k%SN'' 



N''p''dWb + p^WbSN" 



(4.21) 



(4.22a) 



as 



and is given solely in terms of Cauchy data and the undetermined Lagrange multipliers {N}. The condition that the 
lapse contribution must satisfy is more complicated: 



5En 



dSa[N{TlS7T 



(be) 



(ah) 



■[''''1(51^6 + 2fef,(57r["''l - 2v%Sp'' + WSp") 



as 



fcc - per 



TT^'^^^flfc - 2L'^ + 7r[''^lA^, + 27r[''''l66 - 2p''v'l + p''W)6N] . 



(4.22b) 

Note that this is a weak condition, and all determined Lagrange multipliers are properly expressed solely in terms 
of Cauchy data and {N}. 

Just as is the case for GR, one cannot determine the form of these surface terms in general, but must be content to 
determine E^ and Ej^ from (4.22) for each specific situation. As a simple example, one may assume an asymptotically 
flat spacetime in which all of the antisymmetric functions fall off fast enough as r — * cx3 so that they will not contribute 
any surface terms. Assuming the fall-off of the symmetric sector variables as given in [Eol, one finds that there are no 
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additional contributions from (4.22). Evaluating the non- vanishing surface terms: J^j^ dSa N^—g{'-f^'"''>r'^^~ 



assuming the asymptotically Schwarzschild form given in yields the Schwarzschild mass parameter as expected. 

4.5. Completing the Constraint Analysis. It remains to be shown that the Hamiltonian and momentum con- 
straints are preserved under time evolution (they lead to no further constraints), and that they can both be taken 
to be first class, leaving {V'} (4.11) second class. It would then be possible to define Dirac brackets and impose {tp} 
strongly]^ In order to examine the momentum constraints, it is useful to note (and can be explicitly checked) that 
they do indeed generate spatial diffeomorphisms on E: 



{Qi.nain} 



£f[Q]i, 



(4.23) 



and, as this may obviously be extended to any tensor field on S built strictly out of the canonical variables as: 



one finds (by explicit calculation, or using the above results for the last three): 



{V'l[/l]>2[/2]} 

mf],na[n} 



d-'zfif2p— 



zij,£f[f]^0, 



d^z^P2£fm~0, 



- d' 



J2 



(4.24a) 
(4.24b) 
(4.24c) 
(4.24d) 



It has also been checked by explicit calculation that the time evolution of the momentum constraint (4.191) is weakly 
vanishing dtHa ~ 0. 

Thus the only condition required to show that the constraint algebra closes, is that {7i [M] , [iV] } is some linear 
combination of the Hamiltonian and momentum constraints, and therefore vanishes weakly. This is extremely difficult 
to check explicitly since (as noted before) H is not given in a simple form in terms of canonical variables alone, but 
instead depends on Lagrange multipliers which are difficult to solve for in generalp^. This calculation however, need 
not be performed explicitly. 

Considering the dynamics of NOT within the larger arena of hyperspace (the manifold of all spacelike embeddings 
pl|), one may consider the surfaces as embedded in the Riemannian manifold (M, g(" ^-^)) . A particular foliation 
defines a path in hyperspace, the tangent to which is the vector t, as defined in Section 2.1. The hypersurfaces St 



can then be viewed as deformations of an initial hypersurface Sq along the vector field t, and the Hamiltonian must 
then generate the evolution of the fields under these hypersurface deformations. This implies that one must be able 
to write the Hamiltonian in the standard form: d^x [NTL + N'^Ha) |@- It can then be proved from the principle 
of path independence of dynamical evolution , that the Hamiltonian closing relations are constrained to take on 
the form: 



{n[Min[N]} 
{n[M],nb[N'']} 



d'znal^''''\Meb[N] - Neb[M]), 



d^zn£f[M], 



(4.25) 



{Ha [M'^] , Hb [N"] } ^- J^d^zna£f[Mr. 



This principle essentially requires that the data on Eq evolve uniquely to St , regardless of how Sq is deformed into 
Et (i.e., independently of how one foliates the spacetime in-between), or equivalently, that evolution along different 
paths in hyperspace (between identical initial and final points) yield identical results. This also ensures that the 



^This would be accomplished quite simp 
functions of B" and 7"', as determined by (4.11a) and 



V by replac ing u strongly everywhere using (4.11b), and considering both p and as 
Thus (p, u) would no longer be included in phase space, nor would they 



appear anywhere in the Hamiltonian. 

^''Most of these Lagrange multipliers can in fact be ignored in Poisson brackets if Ti is kept in a form very similar to that appearing in 
as their variations will vanish weakly using the algebraic compatibility conditions. The variation of L" (or Ca if ffla has been strongly 
removed from the system) will only vanish if the smearing function is exactly the lapse function A^. This is not the case in general, 
and one would then require the variation of ba with respect to canonical variables in order to compute Poisson brackets containing Ti.; a 
variation that is not easy to compute. 
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collection of all possible Sj's is can be interpreted as describing different slicings of (M, g^^"^^^), and is therefore a 
reflection of the diffeomorphism invariance of spacetimcp]. 



The Hamiltonian for NGT is thus found to be of the form (4.21), consistently generating the field equations and 
compatibility conditions necessary to make Hamilton's equations (4.6) reproduce the field equations (1.3, 1.4, |l.5|) 



derived from the action (1.1) via the Euler-Lagrange equations. The number of degrees of freedom for niNGT may 
now be easily computed using the standard algorithm as given by equation (1.60) of [ p^ . (For now it will be assumed 
that all of the Lagrange multipliers may be solved for, however this issue is discussed further in Section ||.) In this 
case, one finds that the number of configuration space degrees of freedom per spacetime point is: the # of canonical 
coordinates [{Q/} in ( |3.8| )] (13) — the ^ of first class constraints [H and Ha in (4.19)] (4) — the # of second 



class constraints [{-0} in ( 4.11 )] (i x 2) = 8. Since all of the constraints exist in the symmetric sector, two of these 



are the propagating modes of GR, and the remaining six occur in the NGT sector. In the next section, the dynamics 
of pure GR configurations is reduced to that of GR (with a subtlety), and the limit to old NGT (Einstein's Unified 
Field Theory) is discussed. 

5. The reduction to GR and Old NGT 

It is worthwhile at this point to make further contact with GR by considering how one reduces NGT to GR 
configurations, as well as considering the dynamics of old NGT. 

5.1. GR. There is, of course, an important surface in phase space on which all of the antisymmetric components of 
the fundamental tensor are identically zero. If one sets up the system identically on this surface, it should remain 
on it since the dynamics there are identically those of GR (there is a complication, which will be discussed below). 
One could easily just eliminate all skew objects from the action, and derive the ADM results directly. In the context 
of NGT-type theories one would prefer to view the GR dynamics as occurring in the larger phase space of NGT, 
and view the reduction of dynamics as the imposition (by hand) of a constraint on the initial Cauchy surface that 
puts all skew Cauchy data weakly to zero. Consistent dynamics on this surface will require that this constraint be 
preserved in time. 

Considering first the case when a ^ 3/4, one sets up fields on the initial surface such that all skew sector 
Cauchy data ( p°, tt ^"''!, jgb W g) vanish identically; the algebraic compatibility cond itions requi re that v% — and 
^bc — ba — (( A. 2c ) and ( A.2cj ) respectively). The time evolution conditions (( A. 2a ) and ( A.2b| )) (or the equivalent 
in (L8)) guarantee that the antisymmetric components of the fundamental tensor (i3° and respect ively ) will 

remain zero at later times. The skew sector conjugate momenta vanish by the antisymmetric part of (4.1a) and 
( [4.1b|), and remain zero in evolution once the previous results are inserted into ( |4.9| ) through the equations defining 
(j) (4.10). The skew sector is thus consistently eliminated from any future dynamics once all skew Cauchy data is set 
to zero on the initial hypersurface. The symmetric compatibility conditions (A.l) reduce to those of GR, as do the 
symmetric time evolution equations in (4.8) for (what is then unambiguously) the extrinsic curvature of the surface 
k. The Hamiltonian and momentum constraints (4.19a, 4.19bD then take on the familiar forms (i?^^-* = I'^^R^ab)'- 



n 



-2V,[k'^a'Sik]^-2Q5i 



0. 



(5.1a) 
(5.1b) 



Thus given a ^ 3/4, one finds that if at any time, all of the antisymmetric sector canonical variables vanish, then 
the ensuing dynami cs of the system is identically that of GR. The case where a = 3/4 is slightly different, in that 
when one imposes ( [A.2d|) (and the remaining conditions), ba is left undetermined. Then, although one still finds 
that 9t7r[°^l « from ([4.^), the evolution of is determined from dtp"' ~ ^Nir'^'^'^^bh, and is therefore undetermined. 
Similar behavior occurs with the canonical coordinates jab, Wa, as there are nonvanishing contributions from ba to 
(4.9). Thus the skew sector only remains trivial provided that b is chosen to be zero on every surface, and the 
possibility exists that the system may pass through a GR configuration without remaining there, in contrast to the 
case above|^. 



In fact, a close inspection of the arguments in pl| adapted to this case, shows that one only knows (4.25) up to a c-number 
(i.e. independent of the canonical variables). This c-number has been assumed to be zero in this case, as it would not lead to further 
constraints, but to an inconsistent syste m of equations. It is generally believed that the action for NGT generates a consistent set of field 
equations, and conservation laws (A. II) have been derived based on diffeomorphism invariance that imply that this c-number vanishes. 

^2 Note that requiring Stp" ~ StTrl""' 
6a ~ initially. 



on the initial surface is actually a stronger condition than jab^Wa ~ 0, as it implies that 
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5.2. Old NGT. The reduction to old NGT (or Einstein's Unified Field Theory) is accomphshed by setting a ~ 
(setting m — will not affect results, and can easily be left nonzero). The algebraic compatibility conditions in 
this case may be solved for the determined Lagrange multipliers using similar methods to those in ^ , with 

only mild conditions on the components of the fundamental tensor. However (A.3c) is now a constraint related to 
the U{1) invariance that the action will now possess It is straightforward to show that this constraint is first 
class (its Poisson bracket with all other constraints is weakly vanishing) and generates U{1) transformations on S: 
{Wa, e6[p^] [0]} « ea[0]. There is then left five configuration space degrees of freedom per spacetime point, consistent 



with the results for the Unified Field Theory as well as NGT 32 and In the last of these, ad hoc 

constraints are applied in Section 4 so as to remove Wa from the dynamics, guaranteeing a weakly positive-definite 
Hamiltonian for the linearized theory. These constraints cannot in fact be imposed even in the linearized theory jj] , 
since lack of gauge invariance in the full action ^ manifests itself as a non-conserved source for the skew tensor 
field, exciting all five modes even at the linearized level. This is essentially equivalent to the results of Damour, Deser 
and McCarthy Q who avoid discussing the source terms directly by demonstrating that the same effect is generated 
through the coupling of the skew field to a GR background. 



6. A Closer Look at the Compatibility Conditions 



The symmetric sector algebraic constraints are easily solved for in terms of the skew sector Lagrange multipliers 
and Cauchy data. F is determined by (A.4a), Ca is determined by ( A.4b| ), Qa is determined by ( [4.2[), (t° is determined 
by (A.lc), u'\ is determined by (A.le), the trace of which is the constraint ip2- The relation (A. If) can be used to 
uniquely solve f or FL , provided the inverse of 7(°^) exists. In the skew sector, is determined by ( A. 2c ), where 



the trace gives ( |A.3c[ ), determining the Lagrange multiplier W. The remaining equations for A^^ and ba- ( |A.2dD , are 
complicated since one must first replace the symmetric sector Lagrange multipliers, and then solve for the Lagrange 
multipliers in terms of Cauchy data. The behavior relevant to this sector occurs in (A.3cl) alone, leaving the trace of 
( A. 2d ) to determine A^^^. ( A. 3d ) can be written as an algebraic relation for ba solely in terms of Cauchy data (using 
(Ob) and (lA.leD): 



« ebi^I""!] + fcV - fcV + |a7r('"')T^fc + ^[-^"l ^(,,5^ - ia^t"''!/?^!^ - ia^['^''lG[bc]7^'''^W^<i, (6.1a) 
which will be written as (with obvious definitions): 

]50"''6h « (6.1b) 

The issue at hand is whether or not the operator 0°^ is invertible in general. Clearly if the skew sector is weak 
enough and a ^ 3/4, then ^'^"-^^ will dominate the operator, and 0~l will exist. More generally there may exist only 
a subspace on which it is invertible, in which case the component(s) of ba in this space can be solved for, and the 
other(s) are left undetermined. This will also imply that there are constraints corresponding to any undetermined 
components of ba of the form: « 0. One would then need to check that this constraint is preserved in time 
9(5* « 0, in order to determine whether it is possible for the system to remain on this constraint surface. As 
discussed below in the case where the Cauchy data is dominated by the symmetric sector, it appears more likely that 
these are momentary configurations that the system may pass through. 

For the case of mNGT (a = 3/4), one can easily see that the operator 0°"^ is not dominated by the symmetric 
sector, and that in the absence of any antisymmetric components of the fundamental tensor, it disappears altogether. 
Considering (B°- as a perturbation on 7^'^''), all quantities may be expanded to lowest order in powers of 



Pa 

V7, G 



Ba 



[ab) 



= l(ah)B , a 
l(ab), G\^ab] 



= -l[ab] 
-l[ab\ = 



-l(ac)l{bd)l 



cd] 



(6.2) 



If one assumes that ba is of order (B", 7[°''l) (indicated heuristically by ©(skew)), then ( p.lbj ) becomes three constraint 
equations S° w 0. These lead to dtBf^ « for consistency, and one has obtained the dynamics of massive a Kalb- 
Ramond field on a GR background consistent with the linearized treatment of the field equations given in ||] . 

This is not the general case, as one does not have the freedom to assume the behavior of the Lagrange multipliers, 
it must instead be derived from the compatibility conditions. In particular, (3.1b) determines ba in terms of Cauchy 
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data, and keeping the dominant terms, one must find an inverse for the operator O"'': 

Qab ._ ^a^b ^ l^M^^^^j^W ^ ^a^b _ ^{ab)^ . ^ _^ ^a^b ^ (g_3^ 

where in the second form the following have been introduced: 

^lab] ^ ^^abc^^^ 7a:= JVieafccT'"^!, ^ab] ^ ^^e^bcl' , :^ J-e^'^J^.^y (6.4) 

The Levi-Civita tensor density of weight —1: Cabc and +1: e"''"' ~ ^'^"■'^'^ ^'^^'^^ ^'^'^f^ ta&f have been used (e^^^ — eyiz = 
+1), and the notation 7 • 7 = 7°7a has been introduced. (Note that these are not quite duals; the extra numerical 
factor is introduced for convenience.) The inverse of O is found to be: 



0-2 



7 ■ 7 



1 , , B-B-j-j 

^(^f-) ~ + ^'^^b) + 7a7b 

7 • ±) (7 • B)'' 



(6.5) 



(Although the nonperturbative inversion of the operator 0°'' in (|6.1b| ) may be accomplished using similar techniques, 
the general result is not particularly enlightening and will not be given here.) If 7 • 7 = (in which case 7 = Oif one 
assumes that 7^'^'') is nondegenerate) , then assuming that B ^ 0, one finds two constraints coming from projecting 
perpendicular to B, and one relation for ba from that parallel to B. If 7 7^ and either S = or 7-5 = then 
there is one constraint from projecting parallel to 7, and two conditions determining b from the parallel component. 
(The remaining operator when _B — projects out directions transverse to 7). It has not been possible so far to find 
any of these cases which can be maintained in evolution, and so they are expected to be momentary configurations 
and not surfaces on which the system may evolve consistently. 



The solution (3.5) for ba implies that ba ~ O^'E}' is not ©(skew), but in fact C(skew ^). (This does not happen 



in Einstein Unified Field Theory (or old NGT) since a = leaves a term that allows one to solve for ba at lowest 



order.) This in turn shows up in the evolution equations of the antisymmetric sector canonical variables ((4.8) and 



(4.9)), causing them to evolve arbitrarily quickly as the skew initial data is made smaller. It is possible to have 
S sa ©(skew^), in which case b will again reduce to ©(skew), and this behavior is avoided. However, this is a 
condition on the canonical variables alone, and since one could easily set up initial data that does not satisfy this 
condition, it would have to be realized dynamically. At this point it is not known whether this is a reasonable 
expectation, though it seems unlikely that it would occur for any initial data configuration in which the skew sector 
is small^. 

This sort of behavior in the skew sector must be better understood, as the arbitrarily large time derivatives cause 
one to worry about whether GR spacetimes would be unstable in mNGT. It is also not hard to s ee th at the same 
situation may occur near NGT spacetimes as well. As a specific example, consider the solution of ( |6.1a ) specialized 



to spherically symmetric systems. In the general case where both of [B^,^^"^"^^) are non-zero, the solution is: 

Thus one finds that in regions of spacetime where B^ is vanishingly small (for example, perturbative situations, or 
in the asymptotic region of asymptotically flat spacetimes), 61 becomes arbitrarily large, which in turn may cause 
{dtp^,dtWi,dtj2z) to become very large. Although the Wyman sector solution ||, Q (which is becoming the basis 
for most of the phenomenology in mNGT ||3^, |2^, ^) assumes that both B^ and I^i vanish globally, if one 
considers perturbations of these fields on Eq, the above behavior reappears. 

This is essentially the same effect as was found in ||l^ , where the effect of gravitational dynamics on the constraints 
of various derivative coupled vector fields was studied. It was found that constraints on the vector field may be lost 
when GR is considered as evolving concurrently with the vector field (as opposed to the vector field evolving on a 
particular GR background) . This manifested itself as an increase in the number of degrees of freedom in the vector 
field, and singular behavior in the evolution equations when approaching asymptotically flat spacetimes. There one 
finds no evidence of this when considering the vector field dynamics on a fixed GR background. 

This is intimately linked to degenerate solutions of the Euler-Lagrange equations. If one considers linearized 
perturbations to such a solution, the field equations require that the perturbations vanish, whereas the analysis of the 
full dynamics near the degenerate solution may result in quite different behavior. This degeneracy also results in the 
loss of hyperregularity (Chapter 7.4 of ||30[| ) for the Lagrangian system, and equivalence of the Hamiltonian system 

^^Note also that although one might expect that Ti. might impose this constraint on the system due to the presence of terms of the 
form H, it can be shown that these terms identically cancel once a = 3/4 and AjJ^ = A!^;^^ + <5|^6c] is employed. 
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is not guaranteed. However in this case, Hamilton's equations correctly reproduce the Euler-Lagrange equations 
provided the system avoids configurations corresponding to these degenerate solutions. 



Conclusions 

The Hamiltonian formulation of both old (massless) NGT, as well as the massive theory have been given, identifying 
the former as having five antisymmetric sector configuration space degrees of freedom per spacetime point, and the 
latter six. The symmetric sector of phase space has been enlarged over that of GR by a single canonical pair 
representing the fundamental tensor density The existence of two related second class constraints allows one 

to remove these additional coordinates and recover the same number of degrees of freedom as GR in the symmetric 
sector. The reduction of the system to pure GR field configurations has been discussed and the dynamics of GR 
recovered in the case of massless NGT. In the massive theory, nontrivial dynamics of the antisymmetric sector occurs 
unless an undetermined Lagrange multiplier is chosen to vanish identically on every spacelike hypersurface. 

All fields have been decomposed by making use of a surface basis adapted to g^"^^-*, where the surface has been 
assumed to be spacelike (and the time vector timelike) with respect to all of the causal metrics of NGT {gc}. 
Configuration space has been taken to consist of densitised components of the inverse of the fundamental tensor for 
convenience, broken up into the lapse and shift functions and the surface degrees of freedom: {N, N'^ ,'y^°-'''> , B'',7t''''l). 
The lapse and shift functions remain undetermined Lagrange multipliers, and phase space consists of canonical pairs 
that include all of the remaining components of the inverse of the fundamental tensor, as well as an additional 
redundant pair mentioned above. The Hamiltonian has been put into standard form, and the lapse and shift 
functions enforce the NGT Hamiltonian and momentum constraints (Ti. and Ti.a respectively). 

The Hamiltonian for massive NGT has been found to contain determined Lagrange multipliers that become singu- 
lar near particular field configurations. This may occur when the system is very close to any GR field configuration, 
and has been noted to also occur for Cauchy data near the asymptotically well-behaved static spherically symmetric 
solution with nontrivial NGT sector (the massive Wyman solution). This behavior is not uncommon in derivative 
coupled theories, where the rank of the kinetic matrix changes for particular field configurations. At this point in 
time the consequences of this behavior for the massive theory are not well understood. 
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Appendix A. Surface Decomposed Compatibility Conditions and Field Equations 



The surface decomposition of (L4) and ( 1.5d ) is straightforward, however as these will be required weakly, the 
algebraic conditions that determine the Lagrange multipliers, as well as the r ema ining torsion-free condition, have 
been applied to write them in the given form. The definition of Wa given by ( |3.2| ) has also been used throughout. 



A.l. Compatibility Conditions. The symmetric components of (L4) are found to be: 



g.(afc) _ 
-a 



deJ\/^+r-u-2B'^6a, 
-de^ [7] ^"''^ + 7*''''^ (r + u) - 2u"^ 



b 



B"A° 



61 



(A. la) 
(A.lb) 
(A.lc) 
(A.ld) 
(A.le) 
(A.lf) 
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and the skew components: 

e:|^"l = de^[B]" - + B^"^, +7['^''lab - (1 + |a)7('^^)6f, + \a^''''''Wb, (A.2a) 
C^;'^! = -de^[7]['^''l -BV'' + B''c^''+7[''''l(^ + u)-2v"^ (A.2b) 

^[ab] ^ _^(3)j^j[a6] _ ga^b^ ^ B^< + 7t"''l c, + 7^'^'*) A^^ + 7*'^''^ A:?^ + |a(5["5|;l7('^^) (We " 2b^). (A.2d) 

The first two of each of these are dynamical compatibility conditions, and the rest are algebraic conditions that 
determine Lagrange multipliers. It is useful to have the contractions of a few of these: 

Ci^'^) = u - k'^,, (A.3a) 

^j.(afc) ^_^(3)j^j(afc) +^(afc)^^^B6^a _^[ah]^^^^[M^a^^ (^ g^^) 

^J-'^l = vi^) [B]" - aW, (A.3c) 

4"^' = - + B^u^b +7["''lct + (|a - l)7("^)6fc - la^i^^^Wb. (A.3d) 

Note that there is one further algebraic compatibility condition buried in the dynamical compatibility conditions. 
By taking appropriate combinations in order to make use of (2.5) (replacing the variation with a derivative off the 
surface), one can find: 

(tr := r - W^ba + i«7^"''^/9a(Trfc - 26^), (A.4a) 

giving the remaining algebraic condition that determines F. Repeating this last calculation and considering a deriv- 
ative in the surface yields a condition that determines Ca'- 

:= Ca + jafcB" - ia/3aW - \aG^ab]l'''"'\W ^ - 2b,), (A.4b) 

although this is not independent of the algebraic conditions already found. 



A. 2. Algebraic Field Equations and Conservation Laws. The decomposition of the field equations (L5d) 
yields: 

= -de^ M + 3±± - jm^M^^, (A.5a) 

7^(a±) = -de^ [jabB''] + 3(a±) - A^(a±) , (A.5b) 

7^[a±] - ^de^ [W]a + 3[a±] " M^a±] , (A.5c) 

TZab = de^ [K]ab + 3afc - jmH'Iab, (A.5d) 



the second of which has made use of the Jacobi identity ( |2.8| ). For completeness, the components of the mass tensor 
are: 

M^i = F{F - 1) - 2^^'"''>l3apb + ii^7'"'lG[afc], (A.6a) 

M(a±) = + f3c{G[ab]l^''^^ ~ G ^^b)!^''-^) , (A.6b) 

Af(afc) = -{F - l)G(,fc) - iG(,fa)7[^''lG[,d] + 2paPb - 2aaa^ + \i^'''\G,aGbd - GdbGac), (A.6c) 

M[a^] = -Pa + h""'G[fc,]/3a - PciG^ab)!^''^^ - G[afc]7''''), (A.6d) 

^{F- 2)G[,fc] + iG[,b]7[^'^lG[,rf] + 2/3,ab - 2/3^^, - i7'''''(GcfcG,rf - G.aGbd). (A.6e) 
for which a few identities exist: 

Afx± + 7"'A46 = -2(F - 1) - 7'''''G[,d], (A.7a) 

7'"''^Af(^;>)+7'""M[±6] =0. (A.7c) 
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The remaining objects in (A. 5) {3} are: 



3±± 

3(a_L) 
3(ab) 

3[a±] 
3[ab] 



\/i'^[ar +Tu- u%u', + a-{aa - 2c„) + v^v^ - ^a{W)' , 
+ UCa - kabCr'' - (aa - Ca)B''bb + u\{ab - Cfc) + v\bb + v^Xlt, + \aW(Wa - 26a) 



R 



(3) 



(ri\ 

(ab) " '^(a Hfe) + (r + '")fcafc " Oaflb " fccb" a ~ kacU% 
b - jacWl 

fc , „_^„.b , „b L_ „,b 



6„6b - wv'i - jbcv'a + Kd^t - \a{Wa - 2ha)(Wb - 26b) 



\'^^a\W] + \w{aa~Ca) 



-V[^^)[iy],]+2V[:^[6],]+V^^^[A]:, 

+ (r + u)jab + ftcAJjfc + aabb ~ abba ~ icbU^a - ja 



(3)r 



7(3) r 



- jabC + (Cfc - ab)v a + Uba - U J)b - u 



kcbv^a- 



(A.8a) 
(A.8b) 

(A.8c) 
(A.Sd) 

(A.8e) 



One may construct the tensor density (these are, of course, just particular combinations of the field equations): 

Us^'^Gbc + s'^^Gcb) - i(g■"^7^Bc + g^^7^cB - <5^g^^7^cz,), (A.9) 

where Qab = T^ab ~ \%baT^ is the equivalent of the Einstein tensor for NGT. It is a fairly straightforward calculation 
to show that the components Q\_ and Q\ are both algebraic field equations: 



-g7^J 



= \ [\/^3±± + 7"''3ab - im2( V^Afii + 7'^''Ma6)] - ru + uS,%a + JabBV'' + k\x^\ 



■I b^ at 



0' 



(A.lOa) 
(A.lOb) 



-ST^{a±) + B 7?.[ab] 

= V^3(a±) + 3[ab]B'' + Jab [-rs" + ^/^^B^ + 2S''B^6, + 7[''^lac - (1 - 

analogous to the Gauss and Codacci relations respectively in GR. The conservation laws for NGT derived in ||2^ 
written in a general basis are: 

Ba := Ve, [0]^ + i7^BcVe^ [g]^^ = 0, (A.ll) 

guaranteeing the existence of four Bianchi-type identities related to the above algebraic constraints. Note that these 
identities cannot be used here in this form, as the compatibility conditions have been imposed in their derivation. 
Although this will not be necessary here, one could in principle derive the form of (A.ll) that would hold without 
imposing the co mpat ibility conditions, and could therefore be used here to derive the closing relations ( 1.25 ) directly 
from the action ( |l.l| ); using the assumption of diffeomorphism invariance in a slightly more operational manner. 
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